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I. INTRODUCTION 

Bose-Einstein condensates (BECs) of atomic gases 
with spin degrees of freedom have revealed a rich variety 
of static and dynamic properties. The first spinor BEC 
was realized in an optical trap by the MIT group jl| in 
a spin-1 23 Na BEC which was found to have the anti- 
ferromagnetic ground state 0. The same group subse- 
quently demonstrated the formation of spin domains Q 
and quantum tunneling across the spin domains Q . Spin 
exchange dynamics were recently observed in both spin- 
1 and spin-2 0, Q 8T Rb condensates. The thermody- 
namic properties of spinor Bose gas were investigated in 
Ref. 0- 

In the above experiments, the spinor BECs were ob- 
served by destructive absorption imaging after the Stern- 
Gerlach spin decomposition parallel to the magnetic field, 
and therefore the observed quantity was the population 
in each magnetic sublevel. Recently, Higbie et al Q per- 
formed nondestructive imaging sensitive to the direction 
of the magnetization, where the probe axis was taken to 
be perpendicular to the magnetic field. The observed 
quantity was therefore the magnetization perpendicular 
to the magnetic field, and the Larmor precession was ob- 
served. Higbie et al. 9] suggested the possibility of the 
spontaneous magnetization of a spin-1 87 Rb BEC due to 
ferromagnetic interactions. Motivated by that study, in 
the present paper, the spinor dynamics in a ferromagnetic 
BEC are investigated in order to elucidate how sponta- 
neous magnetization of the ferromagnetic spin-1 BEC oc- 
curs in an isolated system in which the total spin angular 
momentum is conserved. 

We consider a spin-1 ferromagnetic BEC at absolute 
zero and assume that all atoms in the BEC are initially 
prepared in the m = magnetic sublevel. Here m de- 
notes the projection of the spin on the quantization axis, 
which we choose to be the z axis. Because of the ferro- 
magnetic interaction, we expect that the magnetization 
grows spontaneously. However, global magnetization is 
prohibited because the total magnetization is conserved 
and should remain zero. We find that the system de- 



velops local magnetic domains of various types, which 
depend on the geometry of the trapping potential. We 
will show that an elongated BEC in the m = magnetic 
sublevel is dynamically unstable against the formation of 
magnetic domains, in which the directions of the spin in 
adjacent domains are opposed. This staggered magnetic 
domain structure is also dynamically unstable, and sub- 
sequently develops into a helical structure due to the fer- 
romagnetic interaction. In a tight pancake-shaped BEC, 
on the other hand, a concentric ring structure of the spin 
is created. 

This paper is organized as follows. Section ITT1 briefly 
reviews the formulation of the spin-1 BEC. Section IIIII 
presents the results of our numerical simulations of the 
time evolution of trapped spinor systems. Section IIVI 
discusses the physical origins of the structure formation 
of spinor BECs found in Sec. IIIII and Sec. [3 provides 
concluding remarks. 



II. FORMULATION OF THE PROBLEM 

The s-wave scattering between two spin-1 bosons is 
characterized by the total spin of the two colliding 
bosons, or 2, and we denote the corresponding scat- 
tering lengths by ao and a<i. The mean- field interaction 
energy can be written as f dr\con 2 + ciF 2 ]/2, where the 
coefficients are given by [lOj 
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with M being the mass of the boson. The system is 
ferromagnetic if c\ < 0, and antiferromagnetic if c\ > 0. 
The density of particles is defined by 
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and the components of the spin vector are given by F x = 
(F + + F_)/2, F y = (F + - F_)/(2i), and 

F, = |^i| 2 -|V-i| 2 (4) 

with 

F + = F*_ = V2(r^o + V'oV'-i), (5) 

where is the mean-field wave function for the mag- 
netic sublevels m = 1,0, and —1. Since c\ is negative for 
the spin-1 87 Rb BEC, its ground state at a zero magnetic 
field is ferromagnetic flOl Hlj. 

In the presence of a magnetic field, the energy of an 
alkali atom is shifted, primarily due to the electron mag- 
netic moment. We take this effect into account up to 
the second order in the magnetic field. The linear and 
quadratic Zeeman energies of the hyperfine spin / = 1 
state are given by p = —ji B B/2 and q = /i|£? 2 /(4Ahf), 
respectively, where ji B is the Bohr magneton and Ahf is 
the hyperfine splitting 12]. In the case of the / = 1 87 Rb 
BEC, we have p/(k B B) ~ -34 pK/G and q/(k B B 2 ) ~ 
3.5 nK/G 2 . The multicomponent Gross-Pitaevskii (GP) 
equations are thus given by 

4r= (-^ +v )* +< «* 

+-^(F + Vl + F_V>-i), (6a) 
dtp ( Ti? \ 

+ci (^^o ± i^±i) • (6b) 

In the following section, we will numerically solve the 
GP equations (|6a|) and (|6b[) in order to examine how the 
spontaneous magnetization occurs under the restriction 
of the conservation of the total magnetization. We will 
take the scattering lengths of 87 Rb, ao = 101. 8ae and 
CL2 = 100. 4<2b, from Ref. 0, where a B is the Bohr radius. 
Following the experiments by Higbie et al. 9], we first 
prepare the ground state ^ g to be in the m = — 1 state 
by using the imaginary-time propagation method Q. 
The spin-polarized condensate is then transferred into 
the m = state as an initial state. In order to simulate 
experimental imperfections in this population transfer, 
we assume very small populations in the m = — 1 com- 
ponent (ip-i(t = 0) = O.Ol^g) in the initial state. Thus, 
the initial state is a state of broken symmetry with very 
small magnetizations in the x and z directions, which 
trigger the spontaneous magnetization. If i/j±i(t = 0) are 
exactly zero, the system develops no spontaneous magne- 
tization within the framework of the mean field approx- 
imation. The time evolution of the system is obtained 
by numerically solving the GP equations JSJ) by the finite 
difference method with the Crank-Nicholson scheme. 

We apply the magnetic field of 54 mG along the trap 
axis, i.e., in the z direction which induces the Larmor 



precession of the spin in the x-y plane at a frequency of 
p/(2irTi) — 38 kHz. Since this frequency is much larger 
than the characteristic frequencies of the system's dy- 
namics, we set p = in the calculation to avoid the rapid 
oscillations. We note that this procedure does not alter 
any of the dynamics of the system other than eliminating 
the Larmor precession of the spin around the z axis. 



III. SPIN DYNAMICS IN TRAPPED SYSTEMS 
A. Spin dynamics in a cigar-shaped trap 

We will first numerically simulate the experimental sit- 
uation of Ref. 9]. Although the trap frequencies in the 
experiment are (u x , uj v , uo z ) = 27r(150, 400, 4) Hz, we as- 
sume an axisymmetric trap with the radial frequency co± 
taken to be equal to the geometric mean co± = ^Juj x uj y 
in order to save the computational task. The number of 
atoms is assumed to be N = 4 x 10 6 . 

Figure [l] shows the time evolution of the m = pop- 
ulation (Fig. HJa)) and that of the column density of 
each spin component (Fig. Gib)). The m = popula- 
tion first decreases due to the spin-exchange interaction 
+ 0— >l + (— 1). It should be noted that the cen- 
tral part of the condensate predominantly converts to 
the m = ±1 state (Fig. Q (b), uo±t = 300) because the 
spin-exchange interaction proceeds faster for higher den- 
sity. This phenomenon was observed in Ref. Q- The 
m = ±1 components then swing back to the m = 
state, and then the dynamical (modulational) instability 
sets in at uo±t ~ 460. The time evolution of the m = 
population in Fig.n( a ) marks the onset of the dynamical 
instability when the large-amplitude oscillation changes 
to a small- amplitude oscillation at uo±t ~ 500. After 
uo ±t ~ 700, irregular spin domains appear, as shown in 
the bottom panel in Fig.Q(b), and each spin component 
as a whole is almost equally populated. Here, the domain 
structure is not static but changes in time in a complex 
manner. We note that the total density n (the top im- 
age in each set in Fig. [2(b)) remains almost unchanged, 
despite the dynamic evolution of each spin component. 
This is due to the fact that the spin-independent inter- 
action energy, con/k B ~ 100 nK, is much larger than the 
spin-dependent interaction energy, c\njk B ~ 1 nK. 

The time evolutions and spatial dependences of the 
three components of the spin vector are shown in Fig. [2 
where the Larmor precession at 38 kHz is eliminated to 
avoid the rapid oscillations in the figure. The time evo- 
lution of each components of the total spin per atom 
F tot /N = J drF/N is shown in Fig. (a). The ini- 
tial state is taken to have a small m = — 1 component 
(^_i = O.Ol^g, i/j g being the m = — 1 ground- state wave 
function), so that the initial values of the spin are given 
by F^/N ~ 0.014, Fy 0t /N = 0, and F^/N ~ -10~ 4 . 
The x component first increases to ~ 0.16 and then de- 
creases. Due to the quadratic Zeeman effect, the x and y 
components of the total spin are not conserved, but the 
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FIG. 1: (a) Time evolution of the relative population of the 
m = component, J dr\i/jo\ 2 / N , and (b) the column density 
distribution of the total density, J dy^2 m ^_ 1 \ip m \ 2 , and those 
of the m = 1, 0, and —1 components, J dy\^m\ 2 , from top to 
bottom. The trap frequencies are oo± = 2ty x 245 Hz and 
uj z — 2tt x 4 Hz. The magnetic field of 54 mG is applied in 
the z direction. The initial state is given by ipi = 0, ipo = 
_4/ g , and ip-i — O.Ol^g, where ip g is the m — —1 
The size of each strip is 600 x 16 in units of 
t 0.69/xm. 
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z component is conserved due to the axisymmetry of the 
Hamiltonian. 

Figure [2 (b) shows each component of the normalized 
spin density F/n, which corresponds to the mean spin 
vector per atom. Only the x component is shown for 
uo ±t < 580 since the y and z components are negligi- 
bly small. Here, we can clearly see that the staggered 
spin domain structure is formed along the z direction at 
co±t ~ 500, due to the dynamical instability. The forma- 
tion of spin domains with staggered magnetization orig- 
inates from the conservation of the total magnetization. 
Since magnetization of the entire condensate in the same 
direction is prohibited due to spin conservation, sponta- 
neous magnetization can occur only locally at the cost 
of the kinetic energy of the domain walls. As mentioned 
above, the projected total spin on the x-y plane is not, 
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FIG. 2: (a) Time evolution of the three components of the 
total spin per atom, F tot /N = J drF/N, and (b) the radial 
and axial distributions of the three components of the mean 
spin vector, F/n. The conditions are the same as those given 
in Fig. Q and the data are taken on the rotating frame of 
reference, in which the 38 kHz Larmor precession in the x-y 
plane is eliminated for the sake of clarity. Each strip in (b) 
indicates the radial (ordinate) and axial (abscissa) distribu- 
tion, where the bottom left corner is the origin (r = z = 0). 
From t — 440 to 580, only the F x /n distributions are shown 
at every oo±At = 10, because the other components are neg- 
ligibly small. After cj±t = 600, F x /n, F y /n, and F z /n are 
shown from top to bottom. The size of each strip is 300 x 6 
in units of (h/mco±) 1 ^ 2 . 



strictly speaking, conserved. However, since the mag- 
netic field is very weak (B = 54 mG), the nonconserving 
spin components remain insignificant (see Fig. [21 (a)). 

The domains are first created around the center of the 
BEC with \z\ < 60/im, since the dynamical instability 
proceeds faster in a higher-density region; afterward, the 
domains gradually extend over the entire condensate. In 
each domain, the length of the mean spin vector is of 
order of unity, and the ferromagnetic state is thus formed 
locally. At uo±_t ~ 600, the y and z components start to 
grow, and the spin domains evolve in a complex manner. 
The fact that the growth occurs first in the x component 
and then in the other components is due to the initial 
condition, in which small magnetization is assumed in 
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FIG. 3: The mean spin vectors F/natr — seen from the —y 
direction, where the vertical axis is the z axis. The conditions 
are the same as those given in Fig. Q The length of the spin 
vector is proportional to \F\/n, and the color represents F x /n 
according to the gauge shown at the top left corner. The spin 
vector is displayed from z = to z = 54/xm in (a)-(c), and 
from z = 35/xm to z = 89/xm in (d). The Larmor precession 
in the x-y plane is eliminated for clarity of presentation. 



the x direction. 

The z dependence of the mean spin vectors F/n at 
r = is shown in Fig. El The staggered magnetic do- 
mains grow from uo±t = 440 to 520 (see Fig. El (a) to (c)). 
The spin vectors in the domains in (c) slightly incline 
in the z direction due to the small initial value of F z . 
Figure El (d) shows a snapshot at an early stage of the 
growth of the y component. The growing y component 
also has a staggered structure, which is shifted from that 
of the x component by one quarter of the wavelength. 
As a consequence, the spin vectors rotate clockwise or 
counterclockwise around the z axis, thereby forming a 
helical structure, as shown in Fig. El (d). Fragmented 
helical structures having typically one or two cycles of 
a helix are also formed in the course of the dynamics 
after uo±t ~ 700, but the lifetime of these structures is 
very short, i.e., < 10/xs. In contrast, the helical struc- 



ture observed in Ref. 9] is not spontaneously formed, 
but "forced" by an external magnetic field. 

The transverse spin distribution, e.g., F x , can be ob- 
served by the nondestructive imaging method in Ref. |?j . 
The staggered domain structure in Fig. El (c) should then 
look like commensurate soliton trains, since the optical 
thickness is high in the red area with a large F x . The 
red and blue areas alternately blink because of the Lar- 
mor precession. The image of the helical structure moves 
along the z axis due to the Larmor precession as observed 
in Ref. Q. 



B. Spin dynamics in a pancake-shaped trap 



Next, we will examine the case of a pancake-shaped 
trap, in which the axial frequency uo z is much larger than 
both the radial frequency uo±_ and the characteristic fre- 
quency of the spin dynamics. In this case, the spin dy- 
namics occur primarily in the x-y plane, in contrast to 
the case of the cigar-shaped trap. We therefore assume 
that the system can be reduced to two dimensions (2D) 
in the x-y plane. The effective strength of the interaction 
in the quasi-2D system is obtained by multiplying Co and 
ci by (STrmuj^h) 1 / 2 [l^. 

Figure 21 shows the snapshots of the total density and 
the x, y, z components of the mean spin vector F/n, 
where co± = cu z /20 = 2tt x 10 Hz, the number of atoms is 
2 x 10 6 , and the strength of the magnetic field is 54 mG. 
The initial state is taken to be = y/l — (0.01) 2 ^ g , 
ip-i = O.Ol^g, and ipi = 0, where ifj g is the ground state 
of the m = — 1 state obtained by the imaginary-time 
method. We then multiply the initial wave function by a 
small anisotropic perturbation 1 + 0.001 Y^)=i ^ e~ r 
in order to simulate experimental imperfection in prepar- 
ing the ground state, where r = (x 2 + y 2 ) 1 ^ 2 and <p is 
the azimuthal angle. Despite this anisotropic perturba- 
tion, the concentric ring structure of the spin (Fig. 01 
uj±t = 26) arises from the dynamical instability, the ori- 
gin of which is the competition between the ferromagnetic 
interaction and the spin conservation. The spatial distri- 
bution of the directions of the spin vector at uo±_t = 26 
is illustrated at the bottom of Fig. 21 The spin vector 
rotates around the z axis at the Larmor precession fre- 
quency, and hence the spin texture displayed as the red 
and blue rings should be observed alternately in time 
by the imaging method in Ref. 9]. The axisymmetry is 
spontaneously broken at uo±_t ~ 50 and the y and z com- 
ponents of the spin vector also grow. After these events 
occur, the axisymmetry is completely lost and compli- 
cated spin dynamics emerge as shown in the fifth row 
(co±t = 100) of Fig. 0] Throughout this time evolution, 
the total density remains almost unchanged, as in the 
case of the cigar-shaped geometry. 
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FIG. 4: The total density and the three components of the 
mean spin vector F/n of the two-dimensional system (from 
left to right panel), where the abscissa and ordinate refer to 
the x and y coordinates in real space. The size of each image 
is 48 x 48 in units of (h/mio±) 1 ^ 2 . The color for the mean spin 
vector refers to the gauge. The bottom panel illustrates the 
direction of the spin at oo±t = 26, where the color represents 
F x /n as shown in Fig. |3| The Larmor precession in the x-y 
plane is eliminated. 



IV. STRUCTURE FORMATION OF SPINOR 
CONDENSATES 

A. Staggered magnetic domains 

In order to gain a better understanding of the mecha- 
nism of the formation of the staggered magnetic domains 
shown in Sec. IIIIl we perform a Bogoliubov analysis of 
our system. For simplicity, we assume a homogeneous 
system with a density of no , and set the wave function as 
V? = e -^^ h (^ + 0) with ji = c n and ^ = (0, y/n^, 0), 



where the three components refer respectively to the am- 
plitudes of the m = 1, 0, and —1 components. Substi- 
tuting this into Eqs. (|6a|) and (|6bj) and keeping only the 
linear terms with respect to 0, we obtain 



dcj) Ti 



ih- 



5±i 



dt 



& ^2 ■ 
"2M V ±P - 



9±1 



+cin (<A±i +<A^i). 
We can solve these equations by expanding <p as 

(f>(r,t) = (u k e i{k - r - Ukt) + WfcC~ i(fc "*'~ & " , * ) ) 



(7a) 



(7b) 



(8) 



Equation (|7a|) has the same form as the equation without 
the magnetic field, and the eigenenergy is given by [HI 

M 



V £ k(£k + 2c n ), 



(9) 



where Sk = h 2 k 2 /(2M). The corresponding Bogoliubov 
eigenvectors and are both proportional to (0, 1,0), 
which describe the density modulation in the m = 
component. It follows from Eq. that if Cq < 0, the 
eigenfrequency uj^ becomes pure imaginary for long wave- 
lengths, and the system collapses. 
Equation (fFEj) has two eigenenergies 



±P+ V (ek + q)(£k + Q + 2cm ), 



(10) 



and the corresponding Bogoliubov eigenvectors take the 
forms 



and 



ul oc (1,0,0), w+oc (0,0,1) (11a) 



a (0,0,1), «- (x (1,0,0). (lib) 



These two modes describe the spin waves that have spin 
angular momenta ±fi, and the energies are shifted by 
the linear Zeeman energies ±p. The quadratic Zeeman 
term only shifts the single-particle energy from Sk to Sk + 
q. From Eq. (|1U|). we find that the eigenfrequencies o;^ 1 
become pure imaginary if g < or 2c\n§ + q < 0. Since 
ci < and g > for the case with the spin-1 87 Rb BEC, 
the dynamical instabilities in the spin wave occur in a 
high-density region. In such a case, the modes become 
most unstable when the imaginary part of becomes 
maximal, i.e., at wave vectors that meet the following 
requirement: 



s k = max(0, -cin - q). 



(12) 



The eigenvectors (|lla|) and (jllbjl have m = ±1 compo- 
nents, and the coherent excitation of these modes arises 
nonzero magnetization in the x-y plane. Even infinites- 
imal initial populations in the m = ±1 components 
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cause exponential growth if the eigenfrequency is com- 
plex, leading to spontaneous magnetization. However, 
if the quadratic Zeeman term is absent (i.e., q = 0), the 
projected angular momentum on the x-y plane [(i^ ot ) 2 + 
(F* ot ) 2 ] 1/2 , as well as F* ot , where F tot = / drF, must be 
conserved, since the linear Zeeman effect merely rotates 
the spin and the other terms in the Hamiltonian are ro- 
tation invariant. In fact, when q = 0, the imaginary part 
of Eq. (|T0|) vanishes for k = 0, which reflects the fact 
that the uniform magnetization of the entire condensate 
is prohibited. Even in this case, local magnetization is 
possible, as long as the total angular momentum is con- 
served, which leads to the structure formation of the spin. 
Thus, the structure formation of the spin reconciles the 
spontaneous magnetization with the spin conservation. 
This is a physical account of why the k ^ modes in 
Eq. (|T0|) become dynamically unstable. 

From Eq. (|12|) . we find that the most unstable wave- 
length in the situation shown in Fig. [2| is given by (see 
also Refs. 



A 



2nh 



2tt 



^-2m(ciriQ + q) 



11/im. 



(13) 



From Fig. (b), we can see that the size of the stag- 
gered domains around z ~ is roughly equal to 13/im, 
which is in reasonable agreement with Eq. (|13|) . The size 
increases with \z\, since the density decreases with \z\. 
The wavelength (|13|) is also consistent with the interval 
between the concentric rings c± 20/im shown in Fig. 2] 

We note that the spin domains observed in Ref. Q 
may also be regarded as a consequence of structure for- 
mation due to spin conservation. In Ref. |3|, the spin-1 
23 Na BEC is prepared at a 1:1 mixture of the m = and 
m = 1 states; then the length of the initial spin vector is 
\/S/2. Due to the antiferromagnetic interaction, the spin 
vector tends to vanish, whereas the z component of the 
total spin must be conserved. (The x and y components 
need not be conserved because of the presence of mag- 
netic field in the z direction.) In the cigar-shaped trap 
of Ref. |3(, the system thus responds to form staggered 
magnetic domains. In fact, the z component of the to- 
tal spin is conserved, whereas the average length of the 
spin vector decreases from \/3/2 to 1/2 by the domain 
formation. 



B. Helical structure 

The spin vector of the helical structure in Fig. 01(d) can 
be written as (e~ iaz , V2,e iaz )/2, where 2n/a is a pitch 
of the helix. The torsion energy per atom is then given 
by ft 2 a 2 /(4m), and the ferromagnetic-interaction energy 
is given by c\n/2. In order to compare the energy of 
the helical structure with that of the staggered domain 
structure, we assume that the spin vector of the latter 
varies in space as (sin V^cos ^,sin 9 ^-)/\/2. The ki- 
netic and ferromagnetic-interaction energies of this state 



become ft 2 a 2 /(8m) and cin/4, respectively. Therefore, 
the helical structure is energetically more favorable than 
the domain structure when |cin| > ft 2 a 2 /(2m). In the 
present case, h 2 a 2 /(2mkB) ~ 0.2 nK and c\n/k& ~ 1 
nK, and hence the change from the domain structure to 
the helical structure is favored. Since the energy is con- 
served in the present situation, the excess energy associ- 
ated with this structure change is considered to be used 
to excite the helical state. Thus the helix in Fig. |3| (d) 
should be regarded as being in an excited state. This ex- 
plains why the helical structure appears only transiently 
in the present situation. If the excess energy can be dis- 
sipated by some means, the helical structure is expected 
to have a longer lifetime. 

In solid state physics, there are two types of magnetic 
domain walls: the Bloch wall and the Neel wall. At the 
Bloch wall spin flip occurs by tracing a helix, while at the 
Neel wall the spin flip occurs in a plane |2fl|. However, 
the domain walls in Fig. [3| (c) are categorized as neither 
type of wall, since the spin vector vanishes in the middle 
of the wall. In the situation given in Fig. [31 the formation 
of the staggered magnetic domains in the x direction is 
followed by the growth of the y component, leading to 
the formation of the helical structure. The helical struc- 
ture may therefore be regarded as the formation of Bloch 
walls. On the other hand, Neel walls may also be formed, 
depending on the initial conditions. 



C. Magnetic field dependence 



The above calculations were carried out for the mag- 
netic field of B = 54 mG, which gives q/\c\n\ ~ 0.01, so 
the effect of the magnetic field was very small. To in- 
vestigate the magnetic field dependence of the structure 
formation of the spin, we performed numerical simula- 
tions for various strengths of the magnetic field. The 
spin domain structure for B = is qualitatively similar 
to that for B = 54 mG. However, the initial, large-scale 
spin exchange shown in Fig. ^ (b) (w±t = 300) is ab- 
sent in the case of B = 0, since the imaginary part of 
half in Eq. (fTU|) for long wavelengths, ~ (q\q + 2cin|) 1 / 2 , 
vanishes when q = 0. The imaginary part for long wave- 
lengths increases with the strength of the magnetic field 
for B < 500 mG. In fact, when B = 200 mG, which 
gives q/\c\n\ ~ 0.1, the initial spin exchange is enhanced 
compared with the case when B = 54 mG. For B = 200 
mG, the size of the spin domain becomes larger, in agree- 
ment with the fact that Eq. ((13)) gives A c± 16 jam. When 
B = 1 G (q/\c\n\ c± 4), the eigenfrequencies cu^r are al- 
ways real, and no spin dynamics occur. This is because 
for large 5, the m = state becomes the ground state in 
the subspace of F^ ot = 0, due to the quadratic Zeeman 
effect. 
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FIG. 5: The three components of the mean spin vector F x /n, 
Fy/n, F z /n (from top to bottom panel in (a) and from left 
to right panel in (b)) for the initial state ipi = a/1 + O.Ol^g, 
ipo — 0, and ip-i = \/l — O.Ol^g, where ip g is the m = — 1 
ground state. The other conditions are the same as those in 
Fig. H for (a) and in Fig. H for (b). 

D. Initial conditions 

We have investigated the structure formation of the 
spin for several other initial states. For example, Fig. [5| 
shows the structure formation obtained for the ini- 
tial state ifji — y/l + O.Ol^g, = 0, and = 
\/l — O.Ol^g, where ijj g is the m = — 1 ground state. 
Spontaneous magnetization from this initial state also 
proceeds with the formation of the staggered domain 
structure for the cigar-shaped geometry, and the con- 
centric ring structure for the pancake-shaped geometry. 
The magnetization occurs in the z direction, since the ini- 
tial spin has a small z component. Thus, we can conclude 
that the spatial spin structure formed in the spontaneous 
magnetization is essentially determined by the geometry 
of the system and does not depend on the initial state. 

V. CONCLUSIONS 

We have studied the spontaneous local magnetization 
of the ferromagnetic spin-1 condensate that is initially 



prepared in the m = state and subject to spin con- 
servation. An initial infinitesimal magnetization of one 
component gives rise to the exponential growth of that 
component due to the dynamical instabilities. As a con- 
sequence, the staggered spin domain structure is formed 
along the trap axis in the case of an elongated cigar- 
shaped trap, where the local mean spin vector under- 
goes the Larmor precession in the presence of a magnetic 
field. The size of the domain structure agrees with that 
obtained by Bogoliubov analysis (Eq. The heli- 

cal structures also appear transiently. In the case of a 
tight pancake-shaped trap, a concentric ring structure is 
formed. 

The underlying physics of the structure formation in 
the magnetization is the competition between the fer- 
romagnetic spin correlation and the conservation of the 
total magnetization. The magnetization of the entire sys- 
tem in the same direction is prohibited. As a consequence 
of the competition, magnetization occurs only locally, re- 
sulting in spatial structures. In solid state physics, in 
contrast, the magnetization of an entire system is possi- 
ble because neither the energy nor the angular momen- 
tum is conserved, due to the interaction of the system 
with its environment. 

In this paper, we have considered only elongated cigar- 
shaped geometry and tight pancake-shaped geometry. In 
other systems, such as the isotropic system and the rotat- 
ing system with vortices, other intriguing spin structures 
or spin textures are expected to emerge in a spontaneous 
manner, because these systems involve dynamical insta- 
bilities of different symmetries. We plan to report the 
results of these cases in a future publication. 
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